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1. Introduction 

Mean-field theory of spin glasses is now almost complete. It took more than thirty years 
from the introduction of a mean-field model by Sherrington and Kirkpatrick (SK) [1] 
before we understood its solution and in particular its physical meaning. The core of 
the mean-field solution was set rather early by Parisi in his replica-symmetry breaking 
(RSB) scheme [2]. Parisi used, however, the replica trick and a formal procedure of 
breaking a symmetry in the (unphysical) replica space when the limit to zero number of 
mathematical replicas is performed. Since then theorists have striven hard to understand 
the physical meaning of the Parisi solution and to find alternative ways of its derivation 
in order to prove its completeness. 

The effort paid off. We now have reached solid understanding of the physics 
behind the RSB solution [3] and there is a mathematical proof of exactness of the 
RSB construction in the SK model [U [5]. This is not the only output of the 
extensive investigation of mean-field spin-glass systems. The statistical methods 
developed by studying the infinite-range spin-glass systems found broad application in 
interdisciplinary fields such as informatics, optimisation and computational complexity, 
econophysics and biophysics and other frustrated complex and open systems. 

Although our understanding of the mean-field theory of spin glasses is in global 
attributes satisfactory, there still remain a few issues that deserve a more detailed and 
specific clarification. One of such questions is the eventual form of the replica-symmetry- 
breaking solution. Any derivation of a stable equilibrium spin-glass state uses the so- 
called discrete RSB scheme with finite-many hierarchical levels of the order parameters. 
Parisi assumed that in the SK model one needs infinite number of RSB hierarchies and 
derived an implicit formula for the free energy of the SK model with a continuous order- 
parameter function. The latter is now considered as the exact solution of the SK model. 
Unfortunately, the Parisi formula for the averaged free energy of the SK model is only 
implicit and no explicit global solution exists. Moreover, the mathematical proof of 
exactness of the RSB construction does not specify whether a discrete or the continuous 
RSB scheme produces the maximal free energy. At least in one model, a Potts spin 
glass, the one-step RSB solution seems to be stable in a finite region of temperatures 
[6]. It is hence important to understand when one should use the discrete RSB scheme 
with a finite number of order parameters and when the continuous limit is appropriate. 

The aim of this paper is to analyse properties of the discrete RSB scheme leading 
to solutions with a finite number of hierarchical levels of the order parameters. The 
emphasis is laid on the way the discrete scheme goes over to the continuous Parisi 
solution in the limit of infinite-many RSB hierarchies. We show that the continuous 
limit is a process during which we reduce the degrees of freedom to a single continuous 
order-parameter function on a compact interval. We derive the continuous limit of the 
stationarity equations maximising the free energy with large but finite numbers of RSB 
hierarchies and thereby we obtain a functional equation for the order-parameter function 
from the Parisi solution. We first derive the equations in the continuous limit generally 
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and then we illustrate the process of building the continuous limit on the asymptotic 
expansion near the de Almeida-Thouless (AT) instability line of the SK model in the 
external magnetic field. We explicitly evaluate the leading asymptotic contribution to 
the Parisi order-parameter function. 

2. Discrete replica-symmetry breaking scheme 

There is no direct way to the Parisi free energy surpassing the concept of discrete 
replicas in one or another way. We can either use the replica trick to handle averaging 
over random configurations of the spin coupling [2J or we can use real replicas with 
which we test thermodynamic homogeneity (independence of boundary and initial 
conditions) of the chosen macroscopic thermodynamic state |7|. In each of these 
approaches one must extend the discrete multiplicity of the replicated phase space to a 
continuous parameter. This can be done only for a specific structure of the replicated 
phase space, namely a vertical tree containing only child replicas within parental ones. 
Communication between replicas of the same generation materialises exclusively via one 
or more antecedent generations of replicas (common ancestors). The phase space forms 
a hierarchical ultrametric structure. It is then natural to start the investigation of the 
RSB construction with the discrete scheme. 

2.1. Hierarchical free energy and order parameters 

A hierarchical character of the phase space of the SK model is expressed in a hierarchy 
of order parameters standing for generations of the replicated spins. Each generation 
is characterised by a pair of numbers mi and xi- The former "geometric" parameter 
expresses a probability with which the original spins are affected by spins from the 
replicated systems and the latter represents the strength with which the original and 
the replicated spins from the Zth hierarchy interact [7] . A phase space with K hierarchies 
(-KTiSB) is characterised in addition to the SK order parameter q also by K pairs {mi, \i) 
with l = l,...,K$ All these parameters are determined from stationarity equations 
for a hierarchical free energy with K hierarchies. 

It is not the set of pairs {mi,xi} that explicitly appears in the hierarchical free 
energy. In fact, we can construct a free energy functional for either xu Am; = m;_i — mi 
or m h Axi = Xi~Xi+i for / = 1, . . . , K with boundary conditions m = 1 and Xk+i = 0. 
Both set of parameters mi and xi form a decreasing sequence. The general formula for 
the hierarchical free energy with the former pairs was introduced in Ref. [8] while the 
latter in Ref. [7j that we use also in this paper. 

The free energy with K hierarchies is characterised by IK + 1 order parameters. 
It is the SK order parameter q and K pairs {mi, Axi}, I = 1, • • ■ ,K. The averaged 

% We use here, in accordance with Ref. [7] , a decreasing sequence of overlap susceptibilities xi instead 
of parameters qi in the Parisi notation. The two sequences are simply related xi = Qk+i-1 — Qo for 
I = l,...,K + l . 
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free-energy density with these order parameters reads 

P K 



1 p ( K \ 

Axk] mi, ... , m K ) = --^ ln2 - - |l - g - A^/J 
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A hierarchical structure of this free energy is evident from the way its interacting part 
lnZ K is constructed. It is the final state in a sequence of partition functions defined 
inductively 

l/mi 



/oo 
Vh z™\ 
-oo 



with the initial condition Zn = cosh 



(2) 



We abbreviated the 



Gaussian differential VX = dX exp{— X 2 /2} /^/2%. 

Free energy ([I]) is a generalisation of 1RSB and 2RSB solutions obtained by 
Parisi. It is a generating functional for all physical quantities of a K -level hierarchical 
solution. The physical values of the order parameters q, Axi, mi, . . . , Axk, m K are 
determined from stationarity equations maximising free energy f K . To represent the 
corresponding stationarity equations we introduce a set of hierarchical density matrices 
in the space of fluctuating random fields A;. We define pi(rj, Xk, ■ ■ ■ , Xi) = Z™ 1 / (i^ mi } A 
where (X(Xi))x l = j^VXi X(X[). We further introduce short-hand notations t 



I anli 



and (t)i(ri; X K ,..., Aj+i) = (pi ■ ■ ■ (pit)x 1 ■ ■ -)x r 
It is now a straightforward task to derive equations for the order parameters from 
the saddle-point equations of functional f K . We obtain 



Q = {{t)l)v > 

Ax* =(((t)li)K) v -(((t)i)K) n , 

mAy _ 4 ((lnZ t _ 1 M,-((lnZ,M, 
' ^ 2(g + Ef =m Ax,)+A X , 



(3) 
(4) 

(5) 



where index I — 1, . . . , K. 

The discrete RSB scheme does not determine a single solution of the original spin 
model, but rather a set of solutions labelled by the number of hierarchies explicitly 
taken into account. Parameter K is hence a free index that is not determined from 
the free energy. Its physical value is fixed by thermodynamic stability. We take so 
many hierarchies of replicas into account till we reach a stable or marginally stable 
solution. Stable solutions with K hierarchies obey K + 1 stability conditions. They are 
a generalisation of the de Almeida-Thouless stability condition of the replica-symmetric 
(K = 0) solution [9]. They reflect non-negativity of eigenvalues of the spin-glass 
susceptibility [10] and read 

2 V 



A i ,(/) = l-/3 2 




i=i 




> o . 



(6a) 
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and 

A i ,(0) = l-/? 2 (((l-t 2 ) 2 )J^>0. (66) 

With the above equations we are equipped with all the necessary tools for finding a 
solution of the SK model. In particular, we can decide whether only a finite number of 
replica generations is sufficient to reach a stable or marginally stable solution or whether 
we must go to infinite-many replicas and the Parisi continuous limit. Parisi deduced 
from 1RSB and 2RSB solutions and confirmed by a truncated model that indeed we need 
infinite number of replica hierarchies. We have recently confirmed this conclusion by 
solving exactly stationarity equations ©-([5]) in the asymptotic region near the critical 
point in zero magnetic field We now extend this asymptotic solution by involving 
the non-zero external magnetic field. 

2.2. Asymptotic solution near de Almeida- Thouless instability line 

Proximity of the instability line naturally introduces a small parameter that we use 
in an expansion of equations (j3J)-((Sj) for the order parameters. If we denote t = 
t&ah\J3(h + r]y/qo)], where q = (tl) v is the SK order parameter in the replica symmetric 
solution, we can define the small parameter to be 

a = (3 2 ((l-tl) 2 ) r) -l>0. (7) 

It measures the distance from the AT line as well as a deviation from the SK solution. 
Since the magnetic field is nonzero, only xi f° r I = 1, ■ ■ ■ , K and Ami for I = 2, . . . , K 
are small. The parameters q and mi are not small unlike the case of zero magnetic field 
We determine the dominant asymptotic behaviour of the corrections to the SK 
solution in the small parameter a for an arbitrary number of hierarchical levels K. 

To derive the leading asymptotic behaviour we must expand Eqs. (j3J) - ([5]) to first 
two nontrivial orders of the small parameter a. We first need to know the two leading 
asymptotic orders of the parameter c[ea = q+J2iLi Axi — 9+Xi to be able to determine 
the leading asymptotic behaviour of Axi and Ami for I > 1- It means, if we want to go 
beyond the one-step RSB solution. 

The asymptotic limit of the RSB solutions with finite numbers of hierarchies 
practically amounts to an expansion in powers of differences Axi- Each integral over the 
random field A; must be expanded at least up to Axf to determine the leading asymptotic 
behaviour of the order parameters. One can rather easily calculate the leading orders 
of Xi an d mi. These two parameters do not depend in the leading asymptotic order on 
the number of hierarchies used and are determined from 1RSB. To see the dependence 
of the order parameters on K means to generate separate equations for single A\i 
with I > 1. We must, however, lift up a degeneracy in the stationarity equations and 
expand them up to Axf- It is a rather tedious task and we accomplished it with the 
aid of the programme MATHEMATICA. The expansion proceeds in the same manner 
we presented in Ref. [11]. We hence do not repeat the detailed steps of the expansion 
but rather summarise the principal findings. 
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T 

Figure 1. The limiting value of parameter m from 1RSB at the AT instability line as 
a function of temperature. The temperature scale was chosen so that T c {h = 0) = 1. 



We must first expand the SK parameter q to the two lowest nontrivial orders in a. 
The solution is then used to determine the lowest asymptotic order of Xi an d We 
obtain mi = m + 0(a) with 



2 (*g(l~ Q% 



and 



a 1 

xi = t^ww^wv, + 0{a2) ' (9) 

These two parameters do not depend on the number of hierarchical levels used. We 
remind that xi = Qea — <?• The values of the temperature and the magnetic field are 
taken from the SK solution. Parameter m is of order unity even at the boundary of 
the spin-glass phase (AT line) where the small parameter a vanishes. Its temperature 
dependence at the AT line is plotted in Figure [TJ 

The other parameter of the 1RSB solution, xii is proportional to the small 
parameter a from Eq. (J7|) and vanishes at the boundary of the spin-glass phase. Its 
ratio Xi/ a a ^ the AT line as a function of temperature is plotted in Figure [2J The ratio 
diverges in zero magnetic field where m = and both xi and m are linearly proportional 
to e = (T c - T)/T c , while a ~ 29 2 p]. 

It is interesting to notice that the asymptotic solution near the AT line reduces to 
1RSB in non-zero magnetic fields. The SK solution can be asymptotically correct in 
the leading order only in zero field where the parameter m\ = 0. In non-zero fields one 
has to go to 1RSB even in the lowest asymptotic order below the instability line. The 
existence of the AT line hence indicates a replica-symmetry breaking. The instability in 
the magnetic field does not specify whether the discrete or the continuous RSB scheme 
applies in the low-temperature phase. Notice also that condition xi > does not 
necessarily indicate a deviation from the SK solution. If mi = then qEA = q+Xi — Isk- 
Only if both parameters Xi an d mi are simultaneously positive the physics of the SK 
solution is changed to 1RSB. 
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Figure 2. Proportionality of the order parameter xi to the small parameter a along 
the AT line. The ratio diverges at the critical temperature T c = 1 as (T c — T) -1 and 
vanishes at zero temperature as T 2 . 



To disclose the leading asymptotic behaviour of each separate parameter Axi and 
Ami for I — 1) • ■ • , K we must go beyond 1RSB and the leading orders in parameters m 
and Xi- It is firstly the fourth order in a in Eq. (JSJ) from which we find that Axi = Xi/K 
and 



ls . , A' • 1 21 



m, = Wi H Am (10a) 

K 

where we added a superscript to specify the number of hierarchical levels used to 
determine the order parameters Xh m i- Further on, we introduced a parameter 
independent of the number of hierarchies Am = m\ — m\. This parameter has an 
explicit asymptotic representation 

/3 2 X-i ((1 - & (2 (1 - 3tlf + 3(4- l)m (8^ + («J - 1) m))) 

Am ~ w^W), (106) 

Both parameters Xi an d Am are linearly proportional to a. The former, however, 
exists already in 1RSB, while the latter first emerges in 2RSB. Since they do not depend 
on the number of hierarchies used and determines a uniform distribution of parameters 
mi for I — 3, . . . , K, we demonstrated that all characteristic features of the asymptotic 
solution near the AT instability line are contained already in 2RSB. What was, however, 
highly nontrivial was to unveil equidistant distributions of both parameters xi an d mi- 
Temperature dependence of the ratio Am /a along the AT line is plotted in Fig. [3j 

Finally we evaluated the instability conditions from Eqs. (l6aj) and (|66l) . They all 
coincide in the leading asymptotic order in a. We find that all the discrete RSB solutions 
are unstable. Instability of the discrete scheme is measured by the small parameter a. 
The SK (K = 0) solution has the instability expressed via the AT condition 

A = l-/5 2 ((l-^) 2 ) r =-a. (11a) 
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Figure 3. The leading-order difference Am factorised by the small parameter a to 
make it of order unity along the AT line. It diverges at the critical temperature T c = 1 
as (T c -T)-\ 




The RSB solutions (K > 1) improve upon stability of the SK solution in that their 
instability is proportional to a 2 . We derived the following explicit expression 

SK Z m + 2 

where Am was defined in Eq. ( 110 fcj) . The difference in the order of magnitude in 
the instability of the SK and RSB solutions is caused by the existence of the non- 
zero parameter m at the AT line where xi = 0. It is hence impossible for the SK 
solution without m to reproduce the exact solution in the non-zero magnetic field even 
asymptotically with a — > 0. 

The stability conditions of the discrete RSB solutions manifest that only the 
continuous limit with K — > oo becomes marginally stable. The instability of 1RSB 
is plotted in Figure HI The leading-order term diverges at the critical temperature 
T c = 1, since the instabilities of the SK as well as of the discrete RSB solutions in zero 
magnetic field are proportional to a oc (T c — T) 2 /T c 2 [TT] . 



3. Continuous replica symmetry breaking 

The above asymptotic solution corroborates the conclusion of earlier calculations in 
zero magnetic field that to reach a stable and consistent solution we need infinite-many 
hierarchical levels in free energy (CQ). It is not generally guaranteed that the limit to 
infinite hierarchies of the discrete RSB scheme must lead to the continuous solution. The 
distribution of parameters Axi = Xi/K calculated in the preceding section explicitly 
manifests that the Parisi solution is the only possible marginally stable solution of the 
SK model in the external magnetic field near the AT line. Moreover, one of us has 
recently proved that the Parisi solution with a continuous order-parameter function can 
always be constructed and is marginally stable in the entire spin-glass phase [12]. Here 
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Figure 4. Leading asymptotic contribution to the instability of 1RSB at the AT line. 
It diverges at the critical temperature as (T c — T)~ 2 . 



we show that the stationarity equations of the solution from Ref. [12] result from the 
continuous limit of stationarity equations (jSJ)-© of the discrete scheme. 

3.1. Homogeneous limit to infinite number of hierarchies 

The basic assumption of the continuous limit of the discrete RSB scheme with the 
number of hierarchies K — > oo is a uniform distribution of the differences A\i- That 
is, they are independent of I and are proportional to K~ x . We then can introduce a 
differential dx = lirnR-_ >00 \\jK. It is not, however, necessary that all parameters Axi 
are equal in the asymptotic limit K — > oo to end up in a continuous theory The 
differences may vary by factors of order unity 

It is convenient to map the hierarchy indices I on interval [0, 1] by introducing a 
continuous variable x = lim^^ 00 (f^ — l)/K. The continuous version of the hierarchical 
free energy results from a process in which we systematically neglect all higher than 
linear orders of the differential dx |13j . It means that in the continuous limit we take 
into account only the second moments of the Gaussian integrations over the auxiliary 
fields A;. We first apply this procedure to the interacting free energy (JTJ and then to 
stationarity equations (jSJ)-©- 

We denote g\ = In Z\. Using the notation from the preceding section we obtain by 
cutting the expansion of Eq. (j2J) at the order O(Axi) 

9l = ln{c .-.)^-^-i- ln{ ^ 1 -. [1 

+^tiA» (glU + | — gt-i + (ft"-, + . (12) 

We denoted g[ = dgi/dh. The derivatives with respect to the magnetic field stand for 
the action of the fluctuating field A; when only the second moment contributes to the 
integral. In the continuous limit we obtain the Parisi differential equation by replacing 
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the hierarchy index I by the continuous variable x 



dg{x, h) x (x) 



dx 



d 2 g(x, h) 
dh 2 



+ m(x) 



dg(x,h,y 
dh 



(13) 



We denoted x(x) = dx(x)/dx. The right-hand side of Eq. fflBl has opposite sign to 
the original Parisi equation, which is caused by a different assignment of the continuous 
variable x to the hierarchy index. Parisi used x p = lim.K->oo l/K. 

To derive the continuous limit of the stationarity equations we must first find a 
reduction of the density matrix pi in the limit of infinite hierarchical levels. Since only 
a linear term in Axi contributes and the density matrix is normalised to unity we have 



pi -> l + mj-iS^AxjJ- . 



(14) 



The operator of the derivative with respect to the magnetic field stands for the random 
variable A;. Integrals over this random variable of functions /(A/, h) weighted with the 
density matrix pi reduce in the continuous limit to 



/ f(\ M\ /en^-L a ' d/((U) A Xl d 2 f(0,h) 
{pif{X h h)) Xi = f(0, h) + mi^Axi — - — + 



(15) 



1 dh '2 dh 2 
We need to evaluate multiple integrals over a number of random fields A;. We 



denote 

fl,i = (Pl+i ■ ■ ■ (Pl+lf)\ l+1 ■ ■ •) 

and using rule ( |T5i) we end up with 

fl,i+l — \Pl+i+lfl,i)\ l+i+1 = fl,i + ^Xl+i+1 



m+isfi+ifli + g/?. 



(16a) 



(166) 



The increment in the second index can again be represented in the continuous limit via 
a differential equation 

dfi(X,h) 



dX 



x(x) 



m{X) —h dh~ + 2 dh 2 



17) 



The solution of the above equation can be represented in form of a "time-ordered" 
exponential with differential operators [12] 

' rX rl 



f x (X, h) = T y exp | J dy X (y) -d- h + m{y)g'(y, h + h)d~ h 



W,h + h) 



h=0 



The ordering operator T y 

rb 



T y exp\ dyO(y)\ = l + J2 d V 



n=l 



yi 



dy2 



Vn-l 



dy n O(yi) . ..0(y n ) 



(19) 



orders products of y-dependent non-commuting operators from left to right in a y- 
decreasing succession. It is a standard tool used in many-body quantum theory to 
represent time-dependent perturbation expansion. It is easy to check that function 
f x (X,h) obeys Eq. (0). 

With the aid of solution (fl8|) we can represent any physical quantity in the 
continuous limit. First among them are the equations for the order parameters. 
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3.2. Stationarity equations and stability conditions 



To derive an equation for the SK order parameter q we simply put f x (0,h) = t(h) 
tanh(/?/i). We obtain from Eq. (CO 



t(X, h) = T y exp | dy x{v) 

Using this representation in Eq. ([3 
q = (t(l,h + r ] ^) 2 



~d\ + m(y)g'{y, h + h)d h 



t(h + h) 



h=0 



(20) 



we reach at a generalisation of the SK relation 



(21) 



To evaluate the right-hand side of Eq. (j3J) we realise that 



Pi (t) 



i-i 



Pi (t)i- 



Axz (t) 



i-i 



from which we find with the aid of integral representation (|T8 

"I 



x( x ) 



T y exp 



dy x(y) 

x(d- h t{x,h v + h)) 2 



d- h + m(y)g'(y, h + h)d h 



h=0/ r) 



(22) 



We denoted h n = h + r\^fq. Equation (!22l) is fulfilled for all x G [0, 1]. It essentially 
determines the functional dependence xi x )- We know from the discrete scheme that 
X{x) > 0. 

The last relation to be rewritten in the continuous limit is equation ([5]). It is not 
difficult to reach a representation 



m(x)(q + x(l) - Xi x )) 



mix) 



, exp 



dy x(y) 



-d\ + m(y)g'(y, h + h)d~ h 



t(x, h v + h)' 



h=0 



(23) 



that now determines the functional dependence m(x). Again from the discrete RSB 
scheme we know that rh(x) < 0. The two functional equations fT22l) and (|23l) allow for 
a trivial solution reducing thereby the RSB quantities to the SK ones. 

The derived equations for the continuous version of the order parameters from 
the discrete RSB scheme enable us to understand how we get rid of one functional 
order parameter. Namely, the function xi x ) appears in all physical quantities only 
under integrals over the index variable x. We hence can redefine the differential 
dx — > dx = dx x(x), since x( x ) > 0. We do not need to know the point-wise 
dependence x( x ) t° determine physical properties of the low-temperature spin-glass 
state. We hence can transform the defining interval x £ [0, 1] to a new one x £ [0, A], 
where X = %(1) < 1. The largest value of x( x ) is the only parameter we need to know 
from this function. It must be determined from a stationarity equation and is related to 
the Edwards-Anderson parameter by an equation q^A = q + X. When we resign on the 
explicit dependence x( x ) we a l so have to disregard Eq. ff22l . Only stationarity equations 
for q and m(x), Eqs. ( 12~TT) and ( 1231) . respectively, remain then relevant. They coincide 
with the equations derived from stationarity conditions imposed upon the Parisi free 
energy in Ref. [12] , 
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Last but not least we have to find the continuous version of stability conditions (16 op 
and (l66j) . It is again straightforward to use the above results and integral 
representation ( fT8l) to arrive at 



A(x)= l-(3 2 {T y exp^ dy X (y) 

dz x(z) m (z)T y exp 
x{d- h t{z, h v + h)) 2 



d\ + m(y)g'(y, h + h)d h 



1 - t(x, h v + hf 



+ 



dy x{y) 

J h=0/ r) 



-dl + m{y)g'{y, h + h)d~ h 



(24) 



It was shown in Ref. [T2] that equality in Eq. f[2~4"|) can be derived from a total derivative 
of Eq. (f23l) with respect to variable x. It means that if equation (!23j) is fulfilled for all 
variables x G [0, 1], both sides of Eq. (I2~4"j) are equal as well. We thereby proved that the 
Parisi continuous RSB solution is marginally stable. The spin-glass susceptibility has 
zero eigenvalue but no negative ones. 



3.3. Asymptotic solution near de Almeida- Thouless instability line 

We now explicitly asymptotically solve the equations for the order parameters of the 
continuous RSB scheme. We know from the discrete version that xi x ) = 1 f° r x £ [0, X\ 
and vanishes elsewhere. We introduce a new dimensionless variable A = x/X that spans 
again interval [0,1]. The physical parameter X serves as an expansion parameter in the 
asymptotic region near the AT instability line. 

In the continuous limit we have two basic equations to solve. It is Eq. ( T2T|) for the 
SK order parameter q and Eq. (1231) from which we determine X and m(A). It can be 
shown that for x < X equation (1221) is a total derivative of Eq. (123]) . There is thus no 
inconsistency if we disregard Eq. ( l22l) as a stationarity equation for the free energy of 
the continuous RSB solution. We know from the preceding subsection that the total 
derivative of Eq. (1231) with respect to x expresses a marginal stability of the continuous 
RSB state. 

The asymptotic solution near the AT line is a polynomial in variable X. We hence 
expand all quantities to a necessary order in this parameter. The T-exponential on the 
right-hand side of Eq. (j23|) must be expanded to X 3 . Further on, the order-parameter 
function becomes also a polynomial in X. The relevant order parameters are then 
expanded as follows 

q =q + Xq[+X 2 q > 2 , (25a) 
m(X) = m + XX m[ . (25b) 

We first use Eq. (121!) to simplify Eq. (1231) and then expand its right-hand side to X 3 . 
Simultaneously we make use of the expansion of the order parameters from Eqs. ( 125 d) 
and (12551) . We obtain a cubic polynomial in A. Coefficients at each power of A must 
vanish and we have three equations for parameters mo,™^ and X. Parameters q[ and 
q' 2 are determined from an expansion of Eq. (|2T|) . The explicit solution was calculated 
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with the aid of the programme MATHEMATICA and reads 

q = (t 2 ) , (26a) 



n 



, 2/? 2 (l-m )(tg(l-tg)) t? 

l-^<(l-^)(l-3^)>,' (266) 

' - ^ (26c) 

Here we denoted 

Q 2 = (l - *g) [(/3 2 (m - 1) (-7t 2 + (51% - 3) m + 5) - mi) tg 

+ (15 (t 2 -l)tl + 2) (3 2 q 12 - 2 (lO*J - 9t 2 + l) (3 2 (mo - 1) <?( . 
The above solution is then used in the equations for X,mo,m' x . These three equations 



read 



(l - tjj) (3 2 (l - 2A/3 2 (m (2t 2 - (37t 4 + 22t\ - l) A/5 2 ) - Am'^ 2 , 
+ (7tl - 3) A/3 2 mgtg) + X/3 2 (l0t 2 + (l05t 4 - 80t 2 + 7)X(3 2 - 2))) , (27a) 
1 - tl) (-211 - 2 (214 - Utl + l) A/5 2 

+m (6 (t 2 - l) Xf3 2 m tl + (-94 + 20t 2 - 3) X{3 2 + 1 - f 2 )))^ , (27b) 

1-t 2 
1 t 



mi(l - 4) + 2/3 2 (l - 3tg) 2 + 3(t 2 - l) /5 2 m (8t 2 + (t 2 - l) m ) ) ) . (27c 



We determine m^ from Eq. fl27cl) , m then from Eq. (12761) . The two parameters we finally 
use in Eq. (I27al) from which we calculate A. To find the leading asymptotic behaviour 
of A we have to expand it in powers of the initial small parameter a measuring the 
depth of penetration into the spin-glass phase. Parameter a was defined in Eq. ((Tj) and 
emerges in Eq. (12 7 aft as an absolute, A-independent term. The solutions for A and mo 
from the continuous RSB scheme then coincide in the leading order in a with the result 
for xi an d m from the discrete 1RSB solution, Eqs. (jUJ) and (jSJ), respectively. The last 
parameter m[ has an explicit representation 

2Am . , 

m[ = - — (28) 

where again parameter Am was already determined within the discrete 2RSB scheme in 
Eq. (11051) . The asymptotic limit of the full solution near the AT line is hence completely 
determined by the parameters from the two-step RSB solution. 



4. Conclusions 



We studied in this paper the behaviour of the replica-symmetry breaking solutions 
in the discrete and continuous schemes. We started with the discrete one with K 
hierarchies and 2K + 1 order parameters q, Axi, mi, ... , Axk, %• Equations for these 
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order parameters are derived from a local maximum of a free energy and were explicitly 
solved in the asymptotic limit to the AT instability line. This calculation served as 
an explicit manifestation of the way the Parisi continuous RSB solution in the non-zero 
magnetic field is approached in the limit K — > oo. We found that unlike in zero magnetic 
field, the SK solution is never, even asymptotically, stable in the non-zero field below 
the AT line. The full solution in the low-temperature spin-glass phase reduces near the 
AT line to the one-step RSB solution. We found that Y$L\ ^Xi = Xi an d mi do not 
depend in the leading asymptotic order on the number of hierarchies K and are exactly 
determined by 1RSB. Further on, we demonstrated that neither Am = K(mi-\ — mi)/ 2 
for I > 2 depends on the number of hierarchies K used. The characteristic parameters 
of the full asymptotic solution are completely set by 2RSB. There is no other parameter 
characterising the asymptotic limit to the AT line. Increasing the number of hierarchies 
in free energy §\§ does not change the values xi — Ya=\ A-Xi an d Am. The new added 
order parameters \u m i for I = 3, . . . , K are equidistantly distributed between the edge 
values < \i < Xi an d m\ — Am < mi < m\ + Am, where Xi an d m\ are calculated 
in 1RSB. These results explicitly prove that the discrete RSB scheme goes over in the 
limit of infinite number of hierarchies to the Parisi continuous RSB solution. 

We analysed the behaviour of the discrete RSB scheme in the limit K — > oo also 
generally. We performed this limit explicitly in the stationarity equations maximising 
the free energy with finite-many hierarchical levels. We derived in this way a set of 
equations for the order parameters in the continuous limit. The equations for the order 
parameters from the discrete RSB scheme goes over in the continuous limit to two 
functional equations for order-parameter functions x(x) and m(x). Since the former 
function comes up only under integrals over the index variable x G [0,1], its point- 
wise behaviour is irrelevant for the physical quantities. The only significant information 
from x( x ) is an integral Jq 1 dxx(x) = X. We hence can disregard the defining equations 
for x(x) and take explicitly into account only equations for the SK parameter q and for 
function m(x). Parameter X is determined from a combination of the two equations. We 
do not lose any information by neglecting the defining equation for x{ x ) that was shown 
to be a total derivative of the equation for m(x). This feature expresses a degeneracy 
of the stationarity equations in the discrete RSB scheme. The equations resulting from 
the continuous limit of the equations from the discrete scheme are identical with those 
derived directly from the Parisi free energy via a saddle point in Ref. [12]. It means 
that the local maximum of the Parisi free energy is a limit of local maxima of discrete 
hierarchical free energies ([I]) when K — > oo. The continuous limit is analytical and all 
physical quantities can be defined and calculated either directly from the Parisi free- 
energy functional of Ref. [12] or from the limit K — > oo of quantities introduced in the 
discrete scheme with free energy (flj) . 

The explicit asymptotic solution of the discrete i^RSB scheme enabled the 
calculation of its thermodynamic stability. We found that in the leading asymptotic 
order of the discrete scheme all the stability conditions ( |6aj) and ( |66j) coincide. Their 
value is negative for any finite number of hierarchies K but approaches zero as K~ 2 . 
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The continuous scheme is then marginally stable with no negative eigenvalue of the 
spin-glass susceptibility. The SK replica-symmetric solution is asymptotically stable in 
the leading order below the critical temperature only in zero magnetic field. In the 
non-zero magnetic field the spin-glass state goes over asymptotically to the one-step 
replica-symmetry-breaking solution (K = 1) that is marginally stable in the leading 
order near the AT instability line. 
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